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Abstract 

We study the percolation behavior of two interdependent scale-free (SF) networks under random 
failure of 1-p fraction of nodes. Our results are based on numerical solutions of analytical expres- 
sions and simulations. We find that as the coupling strength between the two networks q reduces 
from 1 (fully coupled) to (no coupling), there exist two critical coupling strengths q\ and q2, 
which separate three different regions with different behavior of the giant component as a function 
of p. (i) For q > qi, an abrupt collapse transition occurs at p = p c . (ii) For q2 < Q < Qi, the 
giant component has a hybrid transition combined of both, abrupt decrease at a certain p = p J c ump 
followed by a smooth decrease to zero for p < p 3 c ump as p decreases to zero, (iii) For q < q2, 
the giant component has a continuous second-order transition (at p = p c ). We find that (a) for 
A < 3, qi = 1; and for A > 3, q\ decreases with increasing A. (6) In the hybrid transition, at the 
Q2 < Q < Qi region, the mutual giant component jumps discontinuously at p = p 3 c ump to a very 
small but non-zero value, and when reducing p, P^ continuously approaches to at p c = for 
A < 3 and at p c > for A > 3. Thus, the known theoretical p c = for a single network with A ^ 3 
is expected to be valid also for strictly partial interdependent networks. 
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I. INTRODUCTION 



Complex networks appear in almost every aspect of science and technology jl-|l4j. An 
important property of a network is its robustness in terms of node and link failures. The 
robustness of a network is usually characterized by the value of the critical threshold analyzed 
by percolation theory. Recently, motivated by the fact that modern infrastructures are 
significantly coupled together, the robustness of interdependent networks has been studied 
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In interdependent networks, the failure of nodes in one network generally leads to 
failure of dependent nodes in other networks, which in turn may cause further damage to 
the first network, leading to cascading failures and catastrophic consequences. 

The structure of complex networks is frequently non-homogeneous with a broad degree 
distribution. In many cases, the degree distribution obeys a power-law form, and the net- 
works are called scale- free (SF) j^J. Real networks that have been found to be well ap- 
proximated by power-law degree distribution, include between many others, the Internet, 
airline networks, protein regulatory networks, and research collaboration networks (2, 5, 6]. 
Thus, the analysis of interdependent scale-free networks with a power-law degree distribu- 



tion P(k) oc k A is needed. Buldyrev et al. 



developed a framework, based on percolation 



theory, to study the robustness of interdependent networks. Analysis of fully interdepen- 
dent scale-free networks (where all nodes in one network depend on all nodes in the other 
network and vice versa) shows that, the critical threshold is p c > even for A < 3, in 
contrast to a single network where p c = jj]. In general, for fully interdependent networks 
with the same average degree, the broader the degree distribution is (smaller value of A), 
p c is larger 29]. This means that networks with a broader degree distribution become less 
robust compared to networks with a narrower degree distribution. This feature is in contrast 
to the trend known in single non-interacting networks where networks with broader degree 
distribution are more robust. In real world, however, not all nodes in one network depend 
on all nodes in the other network, so it is of interest to study the robustness of two partially 
interdependent scale- free networks. Parshani et al. 30 [generalized the above framework 



29j to study partially interdependent networks. Ref. 30| studied the case of partial cou 



pling where only a fraction q of nodes in each network are interdependent. Their results 
for two interdependent Erdos-Renyi (ER) 3l|, |32] networks show that there exists a critical 



q c , bellow which the system shows a second order percolation transition while above q c a 



first order discontinuous percolation transition occurs. The evolution of such a change from 
first order to second order for SF networks when q changes remained unclear, because the 
behavior of interdependent SF networks is much more complex. 

In this paper, we study the robustness of two partially interdependent scale-free (SF) 
networks under random attack. We assume that only a fraction q of nodes in each network 
are interdependent. We find that for SF networks there are three types of behaviors for 
different q. In addition to first-order transition for large q and second order for small q 
there is a mixed first-second order transition in intermediate q values. Specifically, we find 
(i) As the coupling strength between the two networks, q, reduces from 1 to 0, the giant 
component, of the interdependent networks show three different types of transitions 
with p. For qi < q ^ 1, an abrupt collapse transition occurs. In the range q 2 < q < q±, a 
hybrid transition which is combined of both abrupt and continuous transitions appears. For 
q < q 2 , a continuous second-order transition appears, (ii) The threshold qi which separates 
the first-order and the hybrid transition is equal to 1 for A ^ 3 and decreases with increasing 
A. When q 2 < q < qi, at the steady state of the cascading failures, there exists a p value, 
pi ump , at which the coupled SF networks will suffer a substantial damage due to cascading 
failures but a very small non-zero mutual giant cluster P^ will survive. For p < pi ump , P^ 
will continuously approaches to at p = p c = for A ^ 3 and at p = Pc> for A > 3. 
Thus, the theoretical critical threshold p c = for A ^ 3 for single networks |7| is expected to 
be valid also for strictly partially interdependent networks, (iii) For q < q 2 , the percolation 
transition becomes a regular second order transition, where P^ continuously decreases to 
zero with decreasing p. 



II. CASCADING FAILURES 



Initial failure in one network 



When the system contains interdependent networks, which are several networks fully or 
partially coupled with each other, the initial attack on first network can trigger a systematic 
cascade of failures between the networks (29 1. This can be explained as follows: suppose 
we have a system of two interdependent networks A and B. When, at the initial attack, a 
fraction 1 — p of nodes in network A [A — nodes) are removed since a fraction q of one to one 



bidirectional dependency links exists between A — nodes and B — nodes, so these B — nodes 
which depend on the removed nodes in A are also removed from the network B. Due to 
initial removal, network A may breaks into some connected parts, which are disconnected 
between themselves, called clusters. We assume that only the largest cluster (known as the 
giant component) will function and all the other smaller clusters will become dysfunctional. 
Then the malfunctioning of the nodes in the small clusters of network A will cause failures 
of their counterparts that depend on them in network B, so network B will also breaks into 
clusters, and will cause further fragmentation in network A. This cascade of failures will 
keep going on iteratively until no further failures will occur. 

To theoretically study the pair of coupled SF networks under random failures, we apply 
the framework developed by Parshani et al 3(| to study the cascading failures of partially 
interdependent random networks. Define pa and ps as the fraction of nodes belonging 
to the giant components of networks A and B, respectively. Define ip' n and <p' n as the 
fraction of network A nodes (A — nodes) and network B nodes (B — nodes) remaining, 
and ip n and <fi n the giant components of networks A and B respectively after the cascade of 
failures stage n. Since ip[ stands for the remaining fraction of A — nodes after the initial 
removal, it follows that = p. The remaining functional part of network A therefore 
contains a fraction ipi = ip[pA{ip[)- Because a fraction q of nodes from network B depends 
on nodes from network A, the number of nodes in network B which loses functionality is 
(1 — ipi)q = q[l — ?/4pa(V4)]- Similarly, §' x — 1 — q[l — i/)[pa{i/)[)], among these B — nodes, 
the fraction of nodes in the giant component of network B is <pi = The general 

form of the iterations is 



ip[=p, ip 1 = i)' 1 p A W), 

<P[ = 1- q[l - M)], 01 = fe(^i), 

^ = p[l - q(l -PbM))], ^2 = fe^)-, C 1 ) 
< = p[l - g(l - PJj(&_i))], = <PA«), 



At the end stage of the cascade of failures when nodes failure stops, both networks reach 
a stable state where no further cascading failures happen. According to Eq. it means 
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Ym ~ Ym+li 



(2) 



when m — > oo, since eventually the clusters stop fragmenting. 

Let i)' m be denoted by x and <p' m by y, so we get Vcc = p A (x)x, 0^ = p B (y)y- Applying 
the previous conditions with the last two equations in Eq. (pQ), we obtain the set of equations 



x = p{l -q[l -p B (y)}} 
y = 1 - q[l-pA(x)p]. 



(3) 



Eq. 



30t ] can be solved numerically to get the values of x and y when an analytical 



solution is not possible. This is the case for coupled SF networks, since the generating 
functions of SF network do not have a convergent analytical form, and only an infinite series 
can be obtained. 

Next we introduce the mathematical technique of generating functions for SF networks 



in order to get the analytical forms of p A (x) and p B (x) 
function of the degree distribution is 
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35] . The generating 



G A {z A ) = Y,PA(k)z k A . 



(4) 



Analogously, the generating function of the underlying branching processes is 



H A {z A ) = G' A {z A )/G' A {l). 



(5) 



Random removal of a fraction 1 — p of nodes will change the degree distribution of the 
remaining nodes, so the generating function of the new distribution is equal to t he g enerating 
function of the original distribution with the argument equal to 1 — p(l — z) 33|, l35[. The 
fraction of nodes in A that belongs to the giant component after the removal of 1 — p nodes 
is 



PA ( p ) = l-G A [l-p(l-f A )], (6) 
where f A is a function of p, f A = f A (p), which satisfies the transcendental equation 



f A = G A [l-p(l-f A )]. 
5 



(7) 



For SF networks, the degree distribution is P (k) = ck where A is the broadness of the 



distribution and k m ,i n < k < K. In the case of SF networks [l], 



and 



gam = f: K^r 1 - (^) a - i ]4, (8) 

h—k 



K 

#a(z) = • (9) 



From Eqs. (El)-®, we obtain that 



, _ (l-ZA)(l-C?A(gA)) 

^0° 1-H A (*A) ' / 1Q \ 

, _ (l-z B )(l-G B {z B )) V 7 
^oo - l-H B (z B ) 



where za and satisfy 



SwA = l-q[l-p(l-GA(z A ))], 



1-H b (zb) 

T ^=p[l-qG B (z B )}. 



(11) 



Substituting the generating functions of SF networks into the theoretical frameworks, 
Eqs. (PO)-©, we obtain, using numerical solutions, the theoretical results and compare them 



with results of computer simulations. Fig. 1(a) shows good agreement between the theo 



retical and simulation results for the final giant component ipoo as a function of p for two 
interdependent SF networks under random removal of 1 — p nodes in one network. Three 
cases are studied: (i) A = 2.7, q = 0.95, k min = 2, (k) = 3; (ii) A = 2.7, q = 0.5, k min = 2, 



(k) = 3; and (iii) A = 3.5, q = 0.7, k min = 2, and (k) = 3. Fig. 1(b) shows the cas- 
cading failure dynamics of the giant components left after n cascading stages (denoted by 
ip n ) as a function of number of iterations n, for several random realizations of SF networks 
with A = 2.7, k m i n = 2, (k) = 3 (same parameter values as the numerical calculation), 
N — 1, 280, 000 at p = 0.883 < p c , in comparison with the theoretical prediction of Eq. (CQ). 
Initially the agreement is perfect and when n is getting larger, the random fluctuations in 
topology of different realizations play an important role 
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B. Initial failures in both networks 



When initially a 1 — p fraction of nodes is removed from both networks 37J, |39], the system 
equations d3J) becomes 

x = p{l-q[l-p B (y)p]}, ^ 
y = p{l- q [l _ Pj4 ( x ) p ]}. 

When the degree distribution of the two networks are the same, it follows that Pb(-) = 
Pa(-), x = y and 0oo = ipoo, and the two equations (TT2"]) become a single equation. Further- 
more, using Eqs. ffTOl) and ffTTl) . we obtain 



^ = ^ = ' 1 7» 1 -f» , (13) 
1 — if (z) 



where 2; satisfies 

il^=p{l-g[l-p(l-G( 2( ))]}. (14) 

Eq. (fl4j) is a quadratic equation of g, and only one root has a physical meaning as 
1 = (H(z) -l)[l-q+ y/((l - qY + iggg)] 

p 2(z-l) U ' 1 J 

The maximum of -R(^) corresponds to p c , and 

Pc = max{^( 2c )}' (16) 
where z c is obtained when z — > 1, i.e., ^ = 0, and thus 

maxii?} = lim ^ ~ 1 (1 - g) = H'(l). (17) 

z->l Z — 1 

For two interdependent SF networks, when K — > 00, max{i?} — > 00, so p c = 0. However, in 
the numerical simulations, K can not reach 00, so p c seems greater than 0, but in the theory 
p c = 0. Note that when q — 1, Eq. ffT7|) can yield for max{/?} a finite value since 1 — q = 
and therefore p c can become different from zero as found earlier. 

Now let us relate p c of failure in one network (p°) and both networks failures (p b c ). Our 



previous results 38| show that for two networks (p°) 2 = p b c , so we argue that for two SF 



networks, when p b = 0, it follows that p° = 0. 



III. PERCOLATION BEHAVIOR 



It is known that due to the existence of the interdependence links, when the two-network- 
system is under random attack, the iterative cascade of failures in both networks may result 
in a percolation phase transition that completely fragments both networks when the initial 
fraction of removed nodes is above the critical threshold. When all nodes in both networks 
have 1-on-l dependency links towards their counterpart nodes in the other network (given 
the size of both networks is the same), i.e., q = 1, the percolation phase transition is 



discontinuous and first order 29j; and when the coupling strength q reduces to (which 
becomes the case of a single SF network), a second order percolation transition exists 
However, the change of transition from first to second order for SF networks when q changes 



remained unclear. For coupled Erdos-Renyi (ER) 3l|, |32| networks having Poissonian degree 



distribution a critical point q c exists. For q > q c a first order transition occurs while for 



q < q c a second order continuous phase transition occurs 



30]. 



The percolation behavior of two fully and partially interdependent SF networks, obtained 



from the numerical solutions of Eqs. (EM]), are shown in Fig. 2. Figs. 2(a) and (b) show for 
A = 2.7, the fraction of nodes in the giant component of network A, ^oo, as a function of 
p (fraction of the initially unremoved nodes) for several q values. We can see, as exp ected, 
for SF networks, when q — 1 (fully coupled), the phase transition is first order [29|. This 
means as more and more nodes are initially removed, abruptly, at some value of p = p c , 
the critical threshold, the iterative cascading failure process will completely fragment the 
system. Below p c , there will not exist any cluster of the order of the network size. Thus, 
what still will remain are only very small clusters or single nodes. But just above this critical 
p value, when the failures stops, there exists a giant component in the system. 

When q < 1 but close to 1 (A ^ 3), as p decreases from 1, ip^ first shows a sudden big 
drop similar to q — 1 case, but ip^ does not drop to 0, instead, it drops to a small but 
still non — zero value, which means though the giant cluster in the network suffers a big 



damage, it does not collapse completely (see Fig. 2(b) ). We name the p value where has 



the discontinuous drop to be p> c um P. We mathematically define the p^ Mmp as 

pf mp = {p | max{?M+p) " M~P)}}, (18) 
where + p denotes approaching p from above p, and ~p denotes approaching p from bellow 
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p- 



As p keeps decreasing below pi ump , the small giant component, tp^, smoothly decreases, 



until at p = p c = 0, ^oo wi 
similar to single networks 



l^also reach 0. Thus, the real critical threshold for q < 1 is p c = 
(see the analytical arguments at the end of Section Til Bp . This 



71 



phenomenon can be seen more clearly in Fig. 2(b), which is similar to Fig. 2(a) but the 



y-axis, tp^, is plotted in a logarithmic scale. We see that at p c ump , for q = 0.95, q = 0.9, and 
q = 0.85, the corresponding giant component sizes are reduced from order of 1 by a factor 
in the range of 



(19) 



<M-pf mp ) e [io~ 2 ,io- 4 ], 
^oo(V c ump ) « o(i). 

When p decreases further, ip^ decreases smoothly towards zero for p = (The analytical 
proof is given in Section [HB]). 

This behavior is typical of the behavior of a hybrid-transition, which includes both first 
and second order phase transition properties similar to that found in bootstrap percolation 



40h42| . The giant component first undergoes a sharp jump, which is a characteristic of first- 



order transition, and then smoothly goes to 0, which is a characteristic of a second-order 
phase transition. However, when q is getting smaller, this hybrid-transition phenomenon be- 
comes less apparent, and the percolation phase behavior seems to become, at some threshold 



of q — (72, an ordinary second-order transition. For example, the curve for q = 0.6 in Fig. 2(a 



and Fig. 2(b) seems to suggest a second-order transition since there is no obvious sudden 
drop of the giant component size, instead, it continuously decreases when p decreases. For 
the case of two interdependent ER networks, the system shows either a first order or second 
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3- 39]. 



order phase transition but not a hybrid transition as here 

In network B, which is initially not attacked, similar behavior of the giant component 
0oo can be observed, see Fig. 2(d) However, the difference is that even at p — 0, 0^ does 
not approach to 0, but reaches a finite value. This can be understood due to the partial 
dependency between the networks (q < 1). Even if all nodes in A are removed (p = 0), 
since q < 1, there is a finite fraction, 1 — q, of nodes in B that are not removed and in a SF 

Q 

network any finite fraction of unremoved nodes will yield a giant component [7J. Only in 
the fully coupled (q — 1) case, the mutually connected giant cluster will completely collapse 
at p c > 0. 



A. Estimate of pi ump from Poo as a function of p 

So far we saw (Fig. [2]) that for q 2 < q < qi, as p decreases, the giant component shows 
an abrupt drop similar to a first order transition as Eq. (fTUj) . However, the drop is not to 
Poo = like in a first order transition but to a small finite Poo value. As q decreases, as 
seen in Fig. [3j this drop becomes less sharp and smoother, and tends towards a continuous 
second order transition as in Eq. f )20|) . We analyzed this transition and find that the phase 
transition is like a first order transition with a sharp drop of Poo at pi ump . For q < q 2 , 
the hybrid transition diminishes and the behavior becomes a second order transition with a 
continuous behavior. We are interested to determine the values of q\ and q 2 , which separate 
the three distinct regions. In order to achieve that, we first need to find pi ump . 

To accurately evaluate the values oip J c ump for each q, we compute the number of iterations 



(NOI) in the cascading process which shows a maximum at p> c ump [43(]. The NOI is the 
number of iterative cascade steps it takes the system to reach the equilibrium stage. In the 
simulations, NOI=m is defined by Eq. (J2J), i.e., the step where no further damage occurs. 
But in the numerical solution, ip n is approaching ^/v, only when n — > oo. Here we define 
NOI=m when 

ip m - ip m+1 < f , 

(21) J 

4>m — 4>m+l < £, 

where £ is a very small number. We choose £ = 10~ 16 in this paper, which is equivalent to 
the requirement for the cascading failures to stop in a two-network system when both have 
10 16 nodes. Note that for other very small values of £ the position of p> c ump remains the 
same. 

At the first-order and hybrid order transition point, the NOI has its peak value which 



drops sharply as the distance from the transition is increased 



43| . Thus, plotting the NOI 



as a function of p provides a useful and precise method for identifying the transition point 



pjump Q f hybrid transition. Fig. 2(c) presents such numerical calculation results of NOI. 
The transition point, pi ump , can be easily identified by the sharp peak characterizing the 
transition point. According to the NOI, we define p^ ump as 



pjump = I max {NOI(p)}}. 
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(21) 



From Fig. 2(c) , one can see that the definition of Eq. (1TB]) coincides with the definition 
of Eq. (J2TJ). 



B. Determining q 2 

We know that when the transition is second-order, the order parameter, ij)^, decreases 
continuously. As seen above, ^oo, has a maximum magnitude change at p J c ump , that can used 
to identify the position of p^ um P by Eq. (I2ip . Thus, we can now investigate these maximum 
magnitude changes for different q values at pi ump . In Fig. [3] we plot as a function of 
p only near pl urnp , for different q values ranging between and 1, for several different A 
values. In order to estimate when these changes are discrete and when they are continuous, 
we define 

Woo{ Pi F {q)) 

The rationale for this is as follows. The quantity ip^ ( + pi ump ) is the value of the order 
parameter just before the maximum drop at p > pi ump , and il>oc(~pi ump ) is the value of the 
order parameter right after the maximum drop. In hybrid transition, as we discussed, this 
change in magnitude is large. However, as q becomes smaller, and when finally the transition 
becomes second-order, the change is continuous and the ratio between the magnitudes in 
Eq. f[2"2"j) should becomes 1. Thus, whenever F(q) goes to 0, the corresponding q is q 2 , which is 
the threshold where the hybrid-transition turns into a second-order phase transition. When 
F(q) 0, 

By extrapolating these q positions (where F{q) goes to 0) for different A, we get q 2 as 



a function of A and plot it in Fig. 4(a) Interestingly, q 2 is not monotonic with A but has 
a maximum around A = 2.4. To alternatively identify q 2 , we define the maximum slope as 
function of q as S(q), 

Qf\- r^oo( + P) -1poc(~P)u ■> ft)A s 

S(q) = max{ — \\ P=P ^A- (24) 

When q is below or equal to q 2 , the value of S(q) is very small, representing a continuous 
change of ip^, which is second-order; when q reaches some value, the S(q) has a sudden drop 
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at q2, i.e., the maximum slope becomes dramatically large, representing a sharp change in 
^oo, which is a sign of the occurrence of a hybrid transition. 

By identifying the position of q where the abrupt drop is located, we can also find the 
thresholds q2 which distinguish second-order and hybrid transition. The results shown in 



Fig. 4(c) , Fig. 4(d) , and Fig. 4(e) match very well the results in Fig. 4(a) , supporting our 



method for determining q2 . 



C. Determining q\ 

For coupled SF networks with A ^ 3, only when q = 1 the transition is a first-order, 
which means q\ — 1 for A ^ 3. As A increases above 3, q\ becomes smaller than 1. To 
estimate the q\ values for A > 3, we define according to Eq. ffl9|) the system to have a first 
order transition if satisfies 



^oo(-pf mp ,g)<^. (25) 

Otherwise, it is not a first-order transition. We set here a value a = 10~ n but similar results 
have been obtained for a = 10 -12 and 10~ 13 . We plot qi as a function of A obtained this 



way in Fig. 4(b) 



Now for any given A value, we plot in Fig. [51 p c as a function of q {p c {q))- For A ^ 3, 
only when q — 1 it is a first-order transition, where ip^ abruptly goes to below p c (l); 
When q < 1, it is either hybrid or second-order transition, and t/'oo is strictly only at 
p = for both cases. However, since for the hybrid transition, the giant component becomes 
very small at p } c ump , we can regard this point as an effective p c . For second-order transition, 
although there still exists a p value where there is a maximum change in the magnitude of 
ijjoo, but since ipoo is continuous in all p region, we define p c where ipoo goes to and thus p c 
is always 0. 

For A > 3, the first-order transition happens also for q < 1 and at p c , ipoo jumps to 0. In 
this case p c of the second-order transition and of the hybrid transition is not 0. 
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IV. SUMMARY 



We find that for two SF interdependent networks model with partial dependency q, the 
phase transition behavior of the giant cluster under random attack shows a change from 
first-order (for qx < q < 1) through hybrid transition (52 < q < qi) to a second-order 
phase transition (0 ^ q < q^). in the hybrid transition region, at an effective critical point 
pjump^ giant component V'oo has a sharp drop from finite value to a much smaller, yet a 
non — zero value. The hybrid transition seems to be unique for SF since it does not appear 



in coupled ER networks 



30-|32|. 
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FIG. 1. (Color online) (a) The giant component tpoo as a function of p for coupled SF-SF networks 
system under random removal of 1 — p nodes in one network. SF networks with three different 
parameters are shown (i) A = 2.7, q = 0.95, k m i n = 2, (k) = 3, (ii) A = 2.7, q = 0.5, k m i n = 2, 
(k) = 3, and (iii) A = 3.5, q = 0.7, k m i n = 2, and (k) = 3. The lines represent the theory (Eq. ([3])) 
and symbols are results of simulations, (b) Comparison between theory and simulations of ip n , the 
fraction of the giant component obtained at p = 0.883, which is just below p c , after n stages of the 
cascading failures for several random realizations of coupled SF networks with A = 2.7, k m i n = 2, 
(k) = 3, q = 0.95, and N = 1280000. One can see that for the initial stages the agreement is 
perfect, however for larger n deviations occur due to random fluctuations in the topology between 



different realizations 



36f | . Both simulations and theoretical predictions show a plateau which drops 



to zero, corresponding to a complete fragmentation of the network. Note that some of the random 
realizations converge to a finite mutual giant component, and are not completely fragmented. 
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FIG. 2. (Color online) (a) and (b) Numerical calculations of coupled SF networks with A = 2.7, 
kmin = 2, average degree (k) = (Iva) = (&b) = 3. The size of the giant mutually connected 
component, ^oo, is shown as a function of p for several different values of q. Note that (b) is the 
same as (a) only that the y-axis is in a logarithmic scale. We see that when q = 1, it is a first 
order transition since ipoo goes to zero for p below the jump (pi ump ), but for q = 0.95, q = 0.9, and 
q = 0.85, just below pi ump , ipoo first reaches a small non — zero value, then smoothly goes to zero 
at p = p c = (For analytical proof, see Section [II Bp . This is a typical property of hybrid phase 
transition. For q = 0.6 there seems to be no jump of V'oo and the transition is purely second order, 
(c) The Number-of-Iterations (NOI) 



431 ] to reach the end stage of cascade failure as a function of p. 
(d) Same plot as (a), but for (f)^, which is the giant component of network B, which is not initially 
attacked, (e) and (f) are the same as (a) and (b) but for A = 3.5, and for q = 1, 0.9, 0.8, 0.7 and 
0. 
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FIG. 3. (Color onine) The giant component, ^oo as a function of p for coupled SF networks 
with different values of A, with k m i n = 2, and average degree (k) = 3. Only the critical region 
around the maximum jump of ipoo are shown, for different q values ranging from (most left) to 1 
(most right), with increments of q of 0.01. From these graphs, we can find as the q decreases, ^oo 
becomes more continuous. It is also seen that for large q the sharp jump of starts from small 
but non — zero values to large finite values. This behavior is typical to a hybrid phase transition, 
(a) A = 2.0, the threshold of hybrid transition and second-order transition is q2 — 0.825, so the 
q = 0.82 and q = 0.83 curves are shown with symbols. We can see that the jump in vanishes 
(shown by symbols) when q is reduced, as the phase transition becomes second-order, (b) For 
A = 2.7 and (c) for A = 3.5, the curves in the region of q where the hybrid transition becomes 
second-order are shown by symbols. 
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FIG. 4. (Color online) (a) Values of q% (0) as a function of A for SF networks with average degree 
(k) = 3 and k m in = 2. Note the maximum of q2 at A = 2.4. (b) Values of qi (A) as a function 
of A. Plot of 1/S(q), as a function of q, where S(q) is maximum slope value in the ip vs. p plot, 
are shown for (c) A = 2.0, (d) A = 2.7, and (e) A = 3.5, all SF networks are with average degree 
(k) = 3 and minimum degree k m i n = 2. We can see that the maximum slope values have a sharp 
change at 52 = 0.83, 0.83, and 0.7 for A = 2.0,2.7, and 3.5 respectively, supporting the results in 
(a). 
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FIG. 5. (Color online) The critical threshold p c as a function of q for A = 2.0, 2.7, and 3.5. The 
values of p c are defined as follows: for the first-order transition, p c is where the ipoo jumps to 0; for 
hybrid transition, p , < y ,mp is where the sudden jump of V'oo to a non-zero V'oo occurs; for second-order 
transition, p c is where tpoo goes to 0. For A > 3, we can clearly see three regions of p c . For A < 3, 
qi = 1 and for q < qi sa 0.83, p J c UTnp disappears and p c becomes zero. 
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